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Abstract 

We construct a set of exact ground states with a localized ferromagnetic do- 
main wall and with an extended spiral structure in a deformed flat-band Hubbard 
model in arbitrary dimensions. We show the uniqueness of the ground state for the 
half-filled lowest band in a fixed magnetization subspace. The ground states with 
these structures are degenerate with all-spin-up or all-spin-down states under the 
open boundary condition. We represent a spin one-point function in terms of local 
electron number density, and find the domain wall structure in our model. We show 
the existence of gapless excitations above a domain wall ground state in dimensions 
higher than one. On the other hand, under the periodic boundary condition, the 
ground state is the all-spin-up or all-spin-down state. We show that the spin-wave 
excitation above the all-spin-up or -down state has an energy gap because of the 
anisotropy. 

keywords: ferromagnetic domain wall, spiral state, flat-band Hubbard model, ex- 
act solution, quantum effect, spin-wave, gapless excitation 

1 Introduction 

Domain structures are observed universally in many ferromagnetic systems. If a system 
has a translational symmetry, this symmetry is broken spontaneously by domains. In clas- 
sical spin systems, universal properties of the domain wall have been studied extensively. 
For example in the Ising model on the cubic lattice, Dobrushin proved that a horizontal 
domain wall is stable against the thermal fluctuations at sufficiently low temperatures pQ . 
This structure in the Ising model is also preserved under quantum perturbations. Borgs, 
Chayes and Frohlich proved that the horizontal domain wall on the (/-dimensional hyper 
cubic lattice is stable also against weak quantum perturbations at sufficiently low tem- 
peratures for d > 3 |2j. On the other hand, a diagonal domain wall structure is unstable 
in the Ising model, since some local operators can deform the diagonal domain wall state 
to many other ground states without loss of energy. In the ferromagnetic XXZ model in 
dimensions higher than one, however, no local operator can deform the diagonal domain 
wall ground state to other ground states by the exchange interaction. This fact suggests 
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that the diagonal domain structures are stable at sufficiently low temperatures in suffi- 
ciently high dimensions, even though there is no proof. Alcaraz, Salinas and Wreszinski 
construct a set of ground state with diagonal domain wall structure in the XXZ model 
with a critical boundary field in arbitrary dimensions for an arbitrary spin j^j. Gottstein 
and Werner clarified the structure of ground states in the one-dimensional XXZ model 
with an infinite- volume jlj. Koma and Nachtergaele proved that there is an energy gap 
above any ground states in the one-dimensional XXZ model They also showed an 
interesting result that there exists a gapless excitation above the domain wall ground 
state in the XXZ model in two dimensions jH]. Matsui extended this theorem to the XXZ 
model in arbitrary dimensions higher than one [Zj. Bolina, Contucci, Nachtergaele and 
Starr gave more precise bound for the gapless excitation above the diagonal domain wall 
ground state in the XXZ model jSJ. Bach and Macris evaluate a spin one-point function 
in the domain wall state of the one-dimensional XXZ model by a rigorous perturbation 
method 0. Datta and Kennedy also discussed the existence of a domain wall in one- 
dimensional XXZ models by another rigorous perturbation method ^U]- They show that 
the exchange interaction destroys the domain wall in the antiferromagnetic model, while 
the domain wall exists in the ferromagnetic model at zero temperature. The role of the 
quantum effects should be studied more in many other models. 

Recently, a deformed flat-band Hubbard model with an exact domain wall ground 
state was proposed [TTJ [TJ|. The purpose of this paper is to study excitations above 
the domain wall ground state in this model and clarify whether or not, this model has 
the same spectra as in the XXZ model. The flat-band Hubbard model was proposed 
as a lattice electron model with a ferromagnetic ground state. Some remarkable results 
for ferromagnetic ground states have been obtained in this class of models. Mielke and 
Tasaki have independently shown that the ground state gives saturated ferromagnetism in 
a class of many-electron models on a lattice, which are called flat-band Hubbard models 
|13| ITi] . Nishino, Goda and Kusakabe extended their result to more general models [15 a . 
Tasaki proved also the stability of the saturated ferromagnetism against a perturbation 
which bends the electron band ^H]. Tanaka and Ueda have shown the stability of the 
saturated ferromagnetism in a more complicated two-dimensional model in Mielke's class 
[T7] . Tasaki has studied the energy of the spin- wave excitations in the flat-band Hubbard 
model ^H]- He has shown that the dispersion of the one-magnon excitation is non-singular 
in the flat-band Hubbard model, contrary to the Nagaoka ferromagnetism. The flat- 
band ferromagnetism is believed to be stable against a small perturbation or the change 
of the electron number density [19J. Unlike the ferromagnetic quantum spin model, we 
expect strong quantum effects in the ferromagnetic ground state of the electrons on the 
lattice. The fermion statistics and fully polarized spin configuration imply that this state 
is microscopically entangled with respect to the electron site configuration. Therefore, the 
calculations of the ground state expectation value become more complicated than those 
in the XXZ model in which the ground state can be written in a product state. 

Here, we deform a flat-band Hubbard model by a complex anisotropy parameter q. 
The SU(2) spin rotation symmetry in the original flat-band model is reduced to U(l) 
symmetry in our deformed model. First, we study our model under an open boundary 
condition. The anisotropy \q\ 7^ 1 leads to a localized domain wall with finite width. 
The domain structure is characterized in terms of the local order parameter (S x ), which 
represents the third component of the localized spin at site x. This local order parameter 
takes the same sign within one domain. The domain wall center is a set of sites x defined 
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by zeros of the local order parameter {S x ) = 0. We show the uniqueness of the ground 
state with a fixed magnetization in a half-filled electron number in the lowest energy band. 
We represent (Sl 3 ' ) ) in terms of the local electron density (n x ), and show the profile of the 
ferromagnetic domain wall. We study the low energy excitations in this model. We show 
that there exists a gapless excitation above the domain wall ground state. This excited 
state is constructed by acting a local operator near the domain wall on the ground state. 
We discuss reliability of our results in the infinite- volume limit, although we present our 
result with a finite system size. This property of the domain wall ground state is similar 
to the gapless excitation above the domain wall ground state in the XXZ model as well. 
Next, we study the model under the periodic boundary condition. In this case, either 
all-spin-up or all-spin-down state is allowed as a ground state. We show that a spin-wave 
excitation above the all-spin-up ground state has an energy gap because of the anisotropy. 
This property is similar to the Ising gap in the ferromagnetic XXZ model. 

This paper is organized as follows. In section El we define a deformed flat-band Hub- 
bard model on a decorated (i-dimensional integer lattice. In section El we construct a set 
of ground states and prove the uniqueness of the ground state in a subspace with each 
fixed magnetization. The domain wall structure is shown in terms of the spin one-point 
function. We also obtain a representation for the spin correlation function. In section 
we show the existence of gapless excitation above the domain wall ground state in a suf- 
ficiently large system size. An upper bound on the excitation energy is given in Theorem 
14.11 and Corollary 14.11 In section we consider our model under the periodic boundary 
condition. We estimate an energy gap of the spin-wave excitation above the all-spin-up 
ground state. Finally, we summarize our results in section |U1 



The Hubbard model is a model which represents a many-electron system on an arbitrary 
lattice. In this section, we define a <i-dimensional deformed flat-band Hubbard model 
illustrating its physical meaning. Our model is a generalization of the Tasaki model given 



2.1 Lattice 

The lattice A on which our deformed Hubbard model is defined is decomposed into two 
sublattices 



2 Definition of the Model 



in 



[14j 



A = A U A'. 

A is (i-dimensional integer lattice with linear size L, which is defined 





,2,---,d}. (2) 



A' can be further decomposed to Aj (j = 1, 2, ■ • • , d), i.e. 



d 




(3) 



3=1 



Aj is obtained as a half integer translation of A D to j-th direction, 



Aj := [x + e ij) \x G A } U {x - e {j) \x G A D } , 



(4) 
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where e^' is defined 

e U) :=(0,...,0,i,0,...,0). 

T 

j-th 

We show the lattice in the two-dimensional case in Fig. ^as an example. 



(5) 




Figure 1: Two-dimensional lattice (with L = 3). The white circles are sites in A D and 
the black dots are sites in A'. Electrons at a site can hop to another site if this site is 
connected to the original site with a line or a curve. 



2.2 Electron Operators and the Fock Space 

The creation and annihilation operators for an electron are denoted by cj. a and c XjCr . They 
obey the standard anticommutation relations 

{ c x,ai C y,r} = $x,y$o,T, { c x,ai c j/,t} = = {c x a , C y r }, (6) 

where {A, B} = AB + BA, for sites x, y G A and spin coordinates a, r =|, J,. We define 
no-electron state $ vac by 

for all x G A and a =t, i- We construct a Fock space spanned by a basis 

A,Sca|. (8) 

We also define a number operator n x>a by n x>a = c x a c XjCT whose eigenvalue represents a 
number of electrons at site x with spin a. Note anticommutation relations {cj. CT , cj. CT } = 
z.e. cj. ^cj. a = 0. This relation implies the Pauli principle. We employ the open boundary 
condition, when we consider domain wall ground states. This is realized by c X)CF = if 
\xj\ > L/2 for some j = 1, 2, • • • , d with x = (xi)f =1 . We employ the periodic boundary 
condition, when we consider the spin-wave excitation above the all-spin-up ground state. 
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2.3 Deformed Flat-Band Hubbard Model 

Before we define the Hamiltonian, we introduce new operators a\, a and d xa defined by 



-<f' ° ,/4 E <L#>, + H, - i-' { ° vt E C«v if * e A » 

i=i i=i 
A _1 ct ifxeA' 



and 

A -1 ^ if x G A„ 

-p(a)/A ,\ i if T c A . ' ^ ' 

where q is a complex parameter, A is a positive parameter and p(a) takes +1 if a =| 
and —1 if o —\. And we formally define a xa = and d Xy(J = if \xj\ > L/2 for some 
j = 1,2, ••■ ,d with x = (xi)f =1 . This definitions correspond to the open boundary 
condition for the original electron operators. Note that these a x a and d Xi(T satisfy the 
anticommutation relations, 

We can easily obtain the following inverse relations of © and (fTU|) 

ct _f^ )/4 X>L(^ if^Ao 

c x, CT - ^ i=1 i=i ' ^ 

[A4 i(7 if are A' 

and 

Ag^.o- if a; G A a 

-g- p(CT)/4 4-e(,), CT + \d x , a - q^d x+e{])ta if x e A, ' (13) 

The existence of inverse relations implies that the Fock space is also spanned by another 
basis 



A,BcA . (14) 



ik,t nc • 

VzeA / \x€B / 

This fact is useful to obtain the ground states. 

The definition of our Hubbard Hamiltonian is given by 

H := H hop + H int , (15) 

where H hop and H int defined 

o-=T,l xeA' 

and 

H- mt = U^2n xA n X:l (17) 
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with t,U > 0. The hopping Hamiltonian -f/hop can be written in the following form 



hop 



where 



b x,y 



(42 



[tdW 2 + \q\-V*) 

te -ip(a)e/2 





x,yeA 



if x = y G A 
if x = y G A' 

if x G A', y G A G with [x] < [y] and |x 
if x G A', y G A G with [x] > [y] and \x 
if x, y G A c with [x] > [y] and |x — y| : 
otherwise 



y 



i 

2 



(19) 

with a definition [x] = Ylj=i x j- We parametrize 1 g = \q\e t6 by < 9 < 2ir. Each term 

ti^cj. o-Cy^ in the hopping Hamiltonian represents that an electron with spin a hops from 

site x to site y with a probability proportional to |ti^| 2 . 

Since the interaction Hamiltonian H int represents an on-site repulsive interaction, this 
Hamiltonian can be regarded as a simplification of the Coulomb interaction between two 
electrons. 

Note that this system conserves the number of electron. The total electron number 
operator N e is defined by 

& := £ £ ( 2 °) 

Z£A (7=T,j 

Since the Hamiltonian commutes with this operator, we can set the electron number to 
an arbitrary filling. In the present paper, we only consider that the electron number is 
equal to |A D |, namely, the Hilbert space TC is spanned by the following basis^ 




A,B c A with \A\ + \B\ = |A | 



(21) 



or 




n ~<i ^ 



A,BcA with \A\ + \B\ = \A \ 5> . (22) 

\x€A J \xeB 

Let us discuss the symmetry of the model. First important symmetry is a U(l) 
symmetry. We define spin operators at site x by 



^X^ - = ^ ] C a;,o- 2 Crr ' T ' 

<r,r=U 



(23) 



where (/ = 1,2,3) denote Pauli matrices 



1 

1 



p(2) 



i 



1 
-1 



(24) 



throughout the present paper, we denotes a complex conjugate of a € C by a* and its absolute value 
by \a\. We also denote \v\ to represent a norm of a vector v in rf-dimensional Euclidean space and \A\ to 
represent the cardinality of a set A. 



6 



The Hamiltonian commutes with the third component of total spin operator 

[H, SS] = HS® - S®H = 0, (25) 

with 



7(0 
(3) 

We call the eigenvalue of S^ ( magnetization. We can classify energy eigenstate by the 
magnetization. Note that this U(l) symmetry generated by S^l is enhanced to an SU(2) 
symmetry in the case of q — 1 i.e. Hamiltonian commutes with any component of total 
spin operator. In this case, this model becomes the original flat-band Hubbard model 
given by Tasaki [Til ITSj . Another important symmetry is generated by a product of 
parity and spin rotation defined by 

n = n- 1 = Pexp(m5 t ( 1 t ) ), (27) 

where P is a parity operator defined by Pc x a P = c_ x cr and Pel a P = c_ x a . li transforms 
c XjCr and c xa to C- x ^ and c_ x ^, where a —i if a =| or a =] if a ={. Note the following 

transformation of the total magnetization ILS^II = —S^.. An energy eigenstate with the 
total magnetization M is transformed by n into another eigenstate with the total mag- 
netization — M, which belongs to the same energy eigenvalue. Note that the Hamiltonian 
of the XXZ model with a boundary field h 



d r 



zeAo j=i 



-i 



C(l) c(l) , C(2) o(2) , g + g C(3) o(3) , h c(3) _ u <y(3) 



(28) 



has these two symmetries as well. Our deformation of the hopping Hamiltonian in the flat- 
band Hubbard model is one of the simplest way which preserves these two symmetries. 
If one does not want to add the XXZ Hamiltonian which leads to the ferromagnetism 
trivially, one reaches to our model naturally. 



3 Ground States 

In this section, we obtain ground states of the model with the fixed electron number 
N e = |A C | on the basis of Tasaki's construction method |14j . 

3.1 Construction of Ground States 

The representation of the hopping Hamiltonian in terms of d X(r , 

#hop = t £ £ 4,.4, CT , (29) 

indicates the positive semi-definiteness H hop > 0. The positive semi-definiteness of the 
interaction Hamiltonian H int > is also clear because n x ^ = c xa c x ^ a > 0, then the total 
Hamiltonian is also positive semi-definite 

H = H hop + H int > 0. (30) 
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First, we consider a fully polarized state $^ defined by 



$ T= ( n<t) $ vac. (31) 

We easily verify = from the anticommutativity (jllj) . and therefore $^ is a ground 
state of iJ. Next, we determine all other ground states. 

The conditions that a state $ is a ground state are obviously Hh ov § = and H mt § = 0. 
In other words, 

a Xj(T $ = for all x G A' with a =f, | (32) 

and 

c^Cj^ = for all y G A. (33) 
We expand $ into the following series 



$vac, (34) 



where the summation is taken over all A,BcA with |^4| + \B\ = |A Q |. The first condition 
fl3^|) implies that ip(A, B) does not vanish only for A,Bc A Q . The second condition (J3~3~j) 
for y G A D means ^(A ^0 takes for A n -B 7^ with A,B C A Q . Then we obtain the 
following form: 

* = X>( ff ) ( II S U)*™ (35) 

where the summation is taken over all possible spin configurations a = (cr x )xeA - To 
satisfy the second condition (j22J) for y G Aj (j = 1, 2, • • • , <f), the coefficient holds 

M = ^^-^^"^-^^^(^eO),^)), (36) 

where a x ^ y is spin configuration obtained by the exchange a x and a y in the original con- 
figuration a. This relation implies the uniqueness of the ground state with a fixed total 
magnetization, since two arbitrary spin configurations with same total magnetization can 
be related by successive exchanges of two nearest neighbour spins. Therefore the degen- 
eracy of these ground states is exactly the same as that in the SU(2) symmetric model. 
This degeneracy is also the same as the ground states in the XXZ model pj. 

Note that we can find the "shift operator" S~ which makes the ground state with 
magnetization M from fully polarized ground state $| by acting certain times 

$ M = (S~) Ld - 2M $ T , (37) 
where $m is the ground state with magnetization M. And S~ can be written as 
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3.2 Another Representation of Ground States 

To explore the nature of the ground state, we write the ground state in a more explicit 
way as obtained by Gottstein and Werner in [2J. We define the following electron operator 
creating a superposed state 

a l(0 = Yl %.»4,<7> ( 39 ) 
<t=T4 

where we define a function of x G A and spin a 

„ - A-p(ff)/2 -p(a)[*]/2 

with [x] = X]j=i x j- We define a ground state ^(C) for an arbitrary complex number £ 
by 

*(0 = EC(^)^T = ( II ) $vac (40) 

n=0 \x£A / 

One can see the localization property of the electrons in this representation ^(C)- The 
spin state of an electron at each site is completely determined. Unlike the ground state 
with a fixed total magnetization, one knows the conditional probability of the electron 
spin at a site x. In principle, one can check whether this ground state is realized or not 
by local observations. From this fact, the state \P(C) is expected to be healthy even in 
the infinite- volume limit. We expect that the corresponding ground state to ^(C) is also 
a pure state in the infinite-volume limit, as in the XXZ model. Note that the expectation 
value of an arbitrary local operator in the corresponding ground state ^(C) i n the XXZ 
model is asymptotically equal to that in $a/ for \M\ = 0(1) |S|. We expect that the 
ground state ^(C) i n our Hubbard model has many properties which are the same as 
the domain wall ground state in XXZ model. In our Hubbard model, however, these are 
difficult to be shown, since the state ^(C) defined here in the Hubbard model is not a 
product state unlike the ground state in the XXZ model. 

3.3 Spin One-Point Functions 

Let us now consider expectation values of the spin operators in the ground state ^(C)- 
We denote an expectation value of an operator A in the ground state \I/(C) by (A)^. The 
expectation value of a localized spin at site x is written 

(s u K = i v a, (*(C), 4,^(0) = i v ^ MjCWgg) U1) 
{ x k \fi/ a,T \mo\\ 2 2 CT fl, H*(0II 2 ' 1 ' 

The following anticommutation relations 

{C x ,a, aJ(C)} = ^Vx,aS x ,y, (42) 

for x, y G A D , and 

{C x ,a, <4(C)} = ~Vx,a (q~ 1//L $ x -eW,y + Q 1/4 ^x+e(i) ,y) , ( 4 3) 

for x G Aj and y G A D are useful to calculate the expectation value. These anticommuta- 
tion relations (f4^|l and (|4*H|) yield an equation 

c X! MC) = WMO- ( 44 ) 
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Here, the state ^(C) is defined by 



sgn(x)A JJ aJ(C) ) $ vac if x G A G 



*x(C) = { 



\y^x 



sgn(x - e^)q-^ J] a\{Q + sgn(x + e^)q l ' A J] a\ (C) $ V ac if x G A, 

(45) 

where sgn(x) takes ±1. Then, the expectation value of c* x a c XjT for all x G A in the ground 
state ^(C) can be written in terms of ^ x ((), 

/ f \ * ll V ^'x(C)ll ( AP,\ 

\ C x,ct C x,t/C — Vx^Vx^ \\\Jff£\ || 2 ' V / 

Thus we obtain the representations of spin one-point functions at site x G A in terms of 
an electron number density (n x ){ (n x := n x j + n X) i), 

{S * )c =_ 2 l + |Cg1 2 ' (47) 

m j^c^-ccy (48) 

\ * /c 2i 1 + |C<f| 2 

(5W> C J^k 1 -'^'' 2 (49) 
/C 2 i + |^[x]|2- ^> 

We expect that the electron number density in the ground state ^(C) is almost constant 
on A D or on A' respectively, from the definition of ^(C)- Indeed, in the one-dimensional 
model, we can check this conjecture by the exact bounds [T2] . 

As in the domain wall ground state of the XXZ models discussed in 001201, the two 
domains are distinguished by the sign of the local order parameter (S x )f . The domain 
wall center is defined by zeros of (Sx )c which is located at x with [x] = — logy |C|- The 

function \(n x )^ — \{S X )c\ decays exponentially as x is far away from the center. This 
decay length defines the domain wall width l/log|g|. If the number density is almost 
constant on each sublattice A D or A' as we conjectured, the behaviors of the spin one- 
point functions are not controlled by the number density. In large A limit for real q > 1, 
electrons are completely localized at integer sites, and the spin one-point functions are 
exactly the same as those obtained in the XXZ model defined on A Q . 

For a complex q = \q\e ie , one can see the spiral structure with a pitch angle 9. The 
vector (5 ? x)c := )c)f=i is rotated with the angle [x\9 around the third spin axis 

depending on the site x. Note that this spiral structure of the ground state does not 
exist in the XXZ model, though the complex anisotropy parameter q = e %e is possible in 
the XXZ Hamiltonian. The corresponding model is described in the Tomonaga-Luttinger 
liquid without ferromagnetic order in one dimension. 

The translational symmetry in the infinite-volume limit is broken by the domain wall 
or the spiral structure for finite log | CI - Both symmetries generated by and II are 
broken spontaneously as well. 
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3.4 Spin Correlation Functions 



The spin correlation function can be also represented in terms of the correlation function 
of the local electron number operators 

/q(j)q(l)\ _1 7 lx,a'Pv} r lx,T Vy ^''Pj' \ T 'Vy \t' , > , -r^ 

We can rewrite 

W>>C = (^>C<4°>C7^V- (51) 

if A < 00. If one estimates the correlation function of the local electron number operators, 
one can check the cluster property of the ground state. Actually this can be done for the 
one-dimensional model [T2"] . 



4 Existence of Gapless Excitations 

Here, we show an upper bound of excitation energy in d > 2 for sufficiently large finite 
volume under the open boundary condition. We generalize some parts of Matsui's argu- 
ment for the product ground state in the XXZ model in [7] to those for the non-product 
ground state in the flat-band Hubbard model. We estimate the energy in a trial state 
constructed by acting a local operator on a domain wall ground state ^(C)- 

4.1 Low Energy Excitations 

Here, we show two results for low energy excitations in our model. 

Theorem 4.1 ("Local gapless excitation" above the ground state ^(C) ) I n the d- dimensional 
Hubbard model defined by the Hamiltonian Mty) with d > 2 and the system volume 
|A | = L d , for an arbitrary ( e C with | logy \(\\ < d(L — l)/2 and an arbitrary I with 
1 < I < L + 2| logy |C||/^; there exists a local operator 0\ defined on a compact support 
with a linear size I and a constant F\ > independent of the system size such that 



(0^(0, HOMO) F url (52) 
WOMCW <FlUl ' (5 ) 



and (^(C), Oi^(()) = 0. Moreover, there is an upper bound on the norm of a projected 
state PqOi^>((), where Pq is the projection operator onto the space of ground states. There 
exist constants L\ > 0, R > 1 and F 2 > which are independent of the system size such 
that 

< FS"(L'+ l)R-^ d - forL> L\. (53) 

Here, we describe some physical meanings of Theorem 14.11 We emphasize that the 
excited state in Theorem 14.11 is constructed by acting a local operator which consists of 
finite number of electron operators cj. a and c X)(T . As discussed when we defined the ground 
state ^(C)) one can confirm whether the system takes the ground state ^(C) or not by 
local observations. After one checks the ground state \I/ (C) once, one can obtain a locally 
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deformed state, say Oi^(Q, by a local operation to the system. In this sense, Theorem 14 .11 
claims that one can change the state of the system from the ground state ^(C) by the local 
operation with energy as small as one wants. Particularly, the second result in Theorem 
14.11 guarantees that the deformed state Oi^(() has a non-zero orthogonal component to 
all ground states of the model. This fact implies that the local operation represented in 
Oi is really effective to deform the ground state \P(C) of the system. Therefore, we can 
claim that there exists a gap less excitation above the domain wall ground state ^(C)- 

Remark Theorem 14. II should imply the existence of a gapless excitation in the infinite- 
volume limit under the condition of the fixed electron number, if the corresponding ground 
state in the infinite system to ^(C) were shown to be pure and unique as in the XXZ model 



We can prove the following property of the lowest excitation energy eigenvalue of H 
directly from Theorem 14.11 

Corollary 4.2 (Spectra of a finite system with open boundary condition) Suppose the d- 
dimensional Hubbard model defined by the Hamiltonian \15)) with d > 2. Let E\ be the 
lowest energy eigenvalue of excitation in the model with the system volume |A | = L d 
under the open boundary condition. There exist constants L-2 > and F3 > which are 
independent of the system size such that 



4.2 Expansion in the Original Basis 

Here we prepare for the proof of Theorem 14.21 and 14.11 To evaluate the inner product 
between two states, we represent them in terms of original electron operator cj. a . This 
representation has good properties which help us to estimate expectation values of ob- 
servables. 

4.2.1 Space of Configurations 

To introduce a representation of states in terms of original electron operators, we define 
the decoration of a site x G A by a set 

d 



i=i 

Note A = U x( z\ x, which is not a disjoint union. Also we define the decoration X of a 
subset X C A D by 



To expand the ground state ^(C) in terms of an orthogonal basis, we define the function 
£^(0 for x G A by 




E 1 < F 3 UL 



for L > L 



(54) 






for x G A 
for x G A' 



(55) 
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Since the operator aJ.(C) is written in terms of this function 

4(0 = £ £ ^(04, ff , 

yGx <T=T,J, 

the ground state is represented in 

" " " $vac, (56) 



*(0 = (lIEE^(04J 

\xSA„ j/Sx <r=T,| / 



To represent the ground state in terms of the electron creation operators on sites, we define 
a set Cr of all configurations for the ground state. We define a position configuration for 
the ground state by a one-to-one mapping / : A a — > A with a constraint /(sc) G x for 
each x G A Q . This mapping / selects a site in each decoration x. We denote a set of all 
position configurations for the ground state by Pr. Also we define a spin configuration for 
the ground state by a mapping g : A G — > {til}- We denote a set of all spin configurations 
for the ground state by S. A position configuration / G Pr and a spin configuration 
g G S define a configuration in Cr which is a mapping tp : A c — > A x {f, 1} such that 
</? : x i — > V 9 ( x ) = {fi x )-> 9{ x )) f° r x G A . An arbitrary configuration is also a one-to-one 
mapping. Then the ground state ()56|) is represented as a summation over all configurations 
in orthogonal basis 

*(0 = £ ( n WOi ( ,)Wc. (57) 

Several terms in this summation over all configurations are linearly dependent, and they 
cancel each other. Therefore the summation over all configurations in Cr is reducible. 

4.2.2 Simple Loop 

To obtain an irreducible configuration space, we consider two different configurations ip 
and ip' in Cr. Assume that two terms defined by p and p' are linearly dependent, namely, 
there exists a number C, 



II U-)(Ocl {x) J $vac = C f XI ) ^vac 



(5f 



We say that two configurations <p, p' G Cr are linearly dependent if the terms defined by 
p> and (p' are linearly dependent. This relation implies 

{ip(x)\x G A D } = {(p'(x)\x G A G }, (59) 

otherwise one cannot obtain the relation (|58|) for any number C. Nonetheless, ip ^ ip' 
implies <p{x) ^ (p'{x) for some x G A D . Let us consider a set of sites for two linearly 
dependent ip, <p' G Cr 

X(ip, ip') = {xe A \ip(x) ^ ip'(x)}. (60) 
To study properties of this set, we define several terms. We say that a sequence 

{x m _|_i, X m _|_2, , X m -|_ n } d A Q 

for arbitrary positive even integers m and n is a simple loop with a length n, if x^CiXk+i ^ 
for k = m + 1, • ■ ■ , m + n — 1 and x m+n H x m+ i ^ 0. We say that a configuration 
V 3 = (/jfl 1 ) G Cr for the ground state has a simple loop {x m+ i, x m+ 2, x m +3, • • • ,x m+n ,}, if 
f(x k ) G Xfc+i for A; = m+ 1, m + 2, ■ ■ ■ ,m + n- l and f{x m+n ) G x m +i- 
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4.2.3 Loop Decomposition of Linearly Dependent Configurations 

Here, we show the following lemma. 

Lemma 4.3 Let n be a number of elements of the set X((p, <p') given in f!60|) for two 
linearly dependent configurations <p,<p' G Cr. There exist a positive integer N and some 
positive even integers = m < mi < • ■ ■ < m^ = n for X((p, (p'), such that X((p, <p') can 
be written in a disjoint union of some simple loops 

N 
3=1 

and <p(xk-i) = <p'(xk) f or m j-i + I < k < rrij and <p(x m .) = <p'(x m . 1+ i). Therefore, both 
configurations p and <p' have each simple loop . 

Proof We attach indices to all sites in X(<p, <p') in the following inductive manner. Let x\ 
be an arbitrary site in X(<p, if'). From the relation (J59)) . there exists a site x 2 G X((p, (p') for 
the site x\ such that <p(xi) = <p'(x 2 ). Note X\ 7^ x 2 and xi(lx 2 7^ 0. By the relation (p(x 2 ) 7^ 
(p'{x 2 ) there exists x 3 G X((p,ip') for x 2 such that p(x 2 ) = p\x 3 ) and x 3 ^ Xi,x 2 . There 
exists x 4 G X(ip, if 1 ) for x 3 such that p>{x 3 ) = <p'(x^) and x 4 7^ Xi,x 2 , x 3 . If <p(xi) = ip'(xi), 
then we obtain a simple loop {x%, x 2 , x 3 ,x^} with a length mi = 4 and both configurations 
p and if' have this simple loop. This is because 24 D x\ 7^ 0. If p>{x^) 7^ <p'(xi), we define 
x 5 G X((p,(p') such that (p{x±) = (p f (xs) and 25 7^ xi, x 2 , X3, X4. We assume already 
defined x^-i G X(<p,<p') for an arbitrary natural number k < n with Xfc_i 7^ Xj for any 
i — 1, 2, • • • ,k — 2, such that ^(x?) = ^'(^j+i) for j = 1, — 2. If <p(xk-i) = <p'(xi), 

then we obtain a simple loop {xi,x 2 ,x 3 , ■ ■ ■ ,Xk-i} with a length mi = k — 1 and both 
configurations ip and y?' have this loop. If <p(x^i) 7^ <p'(xi), we define x^ G X((p, cp') such 
that = <p'{xk) and X& 7^ Xi for any i = 1, • • • , k — 1. This can be proved as follows. 

If Xk = Xi for some i — 2, 3, • • • , k — 1, then <£>(x fc _i) = <p'(xk) = <p'(xi) = <£>(xj_i). This 
equality and the definition of a one-to-one mapping imply Xk-i = Xi-i, which contradicts 
the assumption of the inductivity. Also the assumption tp{xk-i) 7^ f'{xi) excludes Xk = 
x\. Thus, Xfc cannot be any element in {x±, x 2 ,x 3 , ■ ■ ■ ,Xk-i}. Note Xk-i fl X& 7^ 0. There 
exists a number m\ < n such that <p(x mi ) = <p'(xi) and x mi flxi 7^ 0. We obtain a simple 
loop {xi,x 2 ,x 3 , ■ ■ ■ ,x mi } with a length m 1; and both configurations p> and have this 
simple loop. If X((p,(p') \ {x 1 ,x 2 ,x 3 , ■■ ■ ,x mi } = 0, then X((p,(f/) = {x 1 ,x 2 ,x 3 , ■■ ■ ,x mi } 
is a simple loop itself with a length n = m\. If X(<p, ip') \ {xi,x 2 , x 3 , ■ ■ • , x mi } 7^ 0, then 
n — m\ > and we continue to attach indices mi + 1, • • ■ , n to the elements in X(<p, <p') \ 
{xi, x 2 , x 3 , ■ ■ ■ , x mi }. By attaching the indices mi + 1, • • • , n, we can continue to identify 
a subset ofX((p, ip') \ {xi, x 2 , x 3 , ■ ■ ■ , x mi } to a simple loop that both configurations have. 
Finally, we can write the set X(tp, ip') as a disjoint union of some simple loops by attaching 
the indices 1,2,3, ••• , n to all the elements in X(ip,ip'). Both configurations tp and ip' 
have each simple loop {x m ._ 1+ i, x mj _ 1+2 , x m ._ l+3 , ■ ■ ■ ,x mj } with a length mj - m^-i, 
and for two sites in each loop we have ip(xk-i) = <p'(xk) f° r m j~i + 1 < k < mj and 
f(x m ) = <p'(x m ._ 1+ i). Thus, we have proved the lemma. I 

4.2.4 Irreducible Configuration Space 

We define the irreducible set of all configurations for the ground state by 

C := {p G Cr I <p has no simple loop}. 
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Here, we obtain the following lemma by proving that the terms given by two configurations 
with a common simple loop are cancelled. 

Lemma 4.4 The ground state is represented in the summation over irreducible set of all 
configurations 

*(o = E Hwoy^nc. (6i) 



Proof From Lemma 14.31 we consider a configuration if G Cr which has a simple loop 
{xi, ■ ■ ■ ,x n } with the length n = 2m. We show that this configuration ip has a counter 
contribution which cancels the contribution from if. For the configuration if, there exists 
a unique configuration ip' with the same simple loop such that X(ip,ip') = {xi, ■ ■ ■ ,x n }. 
In this case, <f(xk) = ip'(xk+i) for any natural number k < n — 1 and f(x n ) = <f'(xi). 
Then, the fermion statistics of the electron operators gives the following relation 



2m 2m 2m 

V(*k)' 

fc=l k=l k=l 



n c Ux k ) - n v(**+i) ~ n 4 



where we define X2m+i = This relation implies C = — 1 in the relation ()58|1. therefore 
the contributions of y> and 99' in ()57|) cancel each other. Now, we rewrite the representation 
(|57|) of the ground state into summation over independent terms. As discussed above, 
all configurations which give linearly dependent terms in the representation (|57|) have at 
least one simple loop and all terms are cancelled. Therefore, the representation (|57|) of 
ground state can be rewritten into a summation over all configurations with no simple 
loop, which consists of only linearly independent terms. I 

The set C of all configurations can be decomposed into a set V of all position configu- 
rations and a set S of all spin configurations, where the set of all position configurations 
is defined by 

V := {/ G I / has no simple loop}, 

as well as C. On the other hand, all spin configurations have no constraint. Therefore, 
after summing over all the spin configurations, the ground state is represented in the 
summation over irreducible set of all position configurations 



*(0 = E 

fev 



II (f/0O.t(C)4(*),T + f/(*y(C)4 (x) J 
.xeA 



$vac (62) 



4.3 Proof of Theorem 14.11 

We define the following hyper cube on the lattice with a positive number I 



Y t := [ < (xi,x 2 , • ■ ■ , Xd) G 



id 



[x] I < , I a;,- 1 < for all j = 1, 2, . . . , d 



z 



y/d\i 



= (1,1,... ,1) HA , (63) 
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where we define z G K. by \C,\ = q z . Note that the sublattice Yi has a linear size I. We 
define a local operator 

^ = \H E ai,^,- (64) 

zeY; <r,T=T,i 

and its deviation from the ground state expectation value 



6S® = 



We define a normalized state ^(C) by 



MOW 2 



5^^(011 



(65) 



(66) 



which is obviously orthogonal to the ground state ^(C)- 
First, we evaluate the norm 



(3) /c(3)\ 



(67) 



Lemma \A. 41 allows us to represent the state in a summation over the configurations. If we 
define an indicator function x by x[t rue ] = 1 an d x [false] = 0, the state can be represented 
in 



y&i feP 



way 



xeA \{y} 



(68) 



where sgn(y) = ±1 is a sign factor coming from the fermion statistics. Note that for an 
arbitrary / 6 V with f(y) ^ y, there exists g E V such that g(y) = y and /|a \M = 
g\x \{y}- Thus we can represent Sy^(C) by 



II «/w,t(04 ( , ))T + (-i) xM e/(-w(c)4 (SB 



$vac + $L- (69) 



The residual state $j_ is orthogonal to the first term in the right hand side as well as any 
ground state \P(C)j since has no term written in a basis of irreducible configurations 
C. The ground state expectation value is represented as 



(sSf>< 



1 \^ \- I£/Qei),t| 2 ~ IO(xiUP TT/lf |2_i_|t I 



(70) 



where P denotes the irreducible set of all the position configurations of the ground state. 



16 



Another important part in the norm of SSy^iC) is 

n4?*(oir = Ei{ E 1+ E 



lf/(*i),T 


2 - 


2 l£/(*a),T 


2 - l£/(x 2 ),ll 2 


I£/(xi),tI 2 + l^/(xi),i 


2 l£/(x 2 ),T 


2 + l£/(x 2 ),ll 2 



> 



E 

/CP I XI GY; 



^/(X!),!! 2 - l£/(xi), 



n(iO(x),Ti 2 +iO(xwi 2 )+ii$ 

xeA 

I^/(xi),t! 2 - K/(xi),|| 2 



X 



±\ 



|2\ 2' 



K/(xi),t| 2 + l£/(xi), 



E 

\xiGY 



IO(xi),tl 2 + ^/(xoal 



X 



n d^Mi 2 + iO(x)ai 2 )- 



XGAo 



Thus, the norm of 5S Yi ^(Q nas the following lower bound 



ii^?»*(oii 2 >e(Ei 



feV UiGY, 



K/(xi),tI 2 - l£/(xi),il 

|C/(xi),t! 2 + l€/(xi),il 



2\ 2' 



+ 



E 

.xi gy 



IO(xi),tI 2 - IO(xi),il 
IO(xi),tI 2 + l£/(xi),il 



-2(4?)C 



n dfe.Ti 2 + i^/(x),ii 2 ) 



XGAo 



> 



EEi 

/GPxiGYi 



K/(xi),t| 2 ~ l£/Ql),il 
K/(xi),t| 2 + IO(xi),|l 



2\ 2 ' 



II (IO(x),tl 2 + IO(x)a 



> 



E 



XI GY; 



| 5 ||[xi]-«|+1/2 _ 


5 


-|[xi]-*|-l/2 


|g||[xi]-«|+l/2 + 


5 


-|[xi]-*|-l/2 



xGA 
\ 2 



EII(i^Ti 2 +iO(xwi 2 )- (7i) 

/GPxGAo 



If we define Gi = (\q\i + \q\ 2) 2 , we obtain 



\\5s^(cw>Gr-t\wo\\ 



d~i\ 



(72) 



Next, we estimate (tfS^tf (C), H5S$H(()). The hopping energy of tffi^tf (£) vanishes 

{6S®*{Q,H^S®*(Q)=0, (73) 

since 5Sy^(() consists of only operators a\ c with i6A„ acting on $ vac . The interaction 
term 

(5S^(0,H iQt 5S^(0) = HtntSgm) = UJ2 Kt^i^OII 2 , (74) 

2/GcY, 



where we define dY t = Y t fl Yf C A'. Note that for any y e <91^, there exist Xi G A \E; 
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and x 2 G Y\ such that x\ D 22 = Thus, we find a bound 

iK,T%i4?*(oir 

\- = y] x[/(gg) = gg] T r / 1 1- ,2, it i2n 

^ (i^),ti 2 + ie/(. 1 wi 2 )(i^),ri a + i^ui 2 ) ii u ^ ),Tl l0( * ul j 

- ,5 (i^),ti 2 + 1 wi 2 xi^ 2 ),ti 2 + U o (lfe ' T ' 2 + lfea|2) 

< ll^(OII 2 r?t . 

- + | g |-M+l+*)2' 

where we have used the indicator function again. For sufficiently large I, there exists 
G 2 > such that 

I < n ld-2 

flolM-i-* + - 2 

Therefore, we have the following inequality 

(6S®*{Q,H6S®*{Q) < G 2 Ul d - 2 \\^(()\\ 2 . (76) 

From (J72j) and (|75jl. we obtain 

(HO,H%(0)<^ur l . (77) 



If we define F\ = G 2 /Gi, we find 

Here, we evaluate norm of Pq^z(C)- First we consider inner product between two 
different ground states 

\mc)MO)\ (7R) 
ii*(coiin*(c)ir 1 ] 

From Lemma 14 A\ we can represent the inner product in the following 

(*(C'),*(C)) _ E^ ec n, eAo WCT WO (79) 



The Schwarz inequality ensures the convergence of the inner product of the normalized 
ground state 

\moMO)\ <1 

ll*(C')llll*(C)ll - ' 

We estimate this inner product with considering the constraints on the configurations. 
We evaluate the product of the norms 

ll*(C')il 2 il*(C)U 2 

= E n(i4(,),Ti 2 + i4(,ui 2 )(ie/ 2 (,) > Tl 2 + ie/ 2( ,ul 2 )- (so) 

fij2ePxeA 
Here, we abbreviate ( and (' by 

£/(x),<t = £f{x),a{(), £f{x),cr = £/(x),<r(C )• 
18 



We evaluate each term in this summation. A term with arbitrary fi, f 2 G V has a lower 
bound 



v i 

^ 2 

h,hev 



II (l^ ( ,),tl 2 + l^(,),il 2 )(l^(.),Tl 2 + \^n\ 2 ) + ih - / 2 ) 



fuhev V ^eA 

2. 



X 



h,hevxeh„ 

ISiW>t ^/i(*)-t + ^/i(«U ^A(*ull^/2(*)A(*).T* + ^/ 2 (x),^/2(x),4*l> 
where a function R(y 1 ,y 2 ) is defined for arbitrary yi,y2 G A by 



(IC, T l 2 + ICJ 2 )(IWI 2 + IWI 2 )x(yi-2/2) 



(81) 



In the practical calculation, we can check R(yi,y 2 ) > 1 for arbitrary ?/i,?/2 £ A. Also 
the Schwarz inequality for the linearly independent two vectors (Cy,h^y,l) anc ^ (^'y,V^y,0 
ensures this relation. We define a function 

R(x) = min R(f 1 (x)J 2 {x)), 

which is also larger than 1. Therefore, we have an upper bound of the inner product 
between the normalized ground states 



*(oiiii*(oii - 11 { } • 



(82) 



xeAo 



Each factor R(x) 1 with a fixed [x] is a constant less than 1, since the both functions 
£,f(x)AO and ZfwAO depend on x only through [x] = Yl=i x j for an Y f € V. We 
define R = ma.x xe \ o R(x), then we find 



iwco^coi^^-^ii^coiiii^oii- 



(83) 



Next, we evaluate an inner product between 5Sy^(() and another ground state. First, 
we evaluate 

z /ePlieyi ?/(xi),t ?/(*i).T + ?/(xi),i xGAo 
Then, we have 



m%s^(o)\ 4E 1 



e n (^)/^)-t + ^w,i*^):i) 

/e-p xeA 



~-G 3 i d \mOM0)\, 



(85) 
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where we define G 3 by YuxieY ^ = G?,l d - Also we obtain another upper bound 

(m,S?M0)<lG 3 l d \MC)\\ 2 . (86) 
Therefore, the inner product between SSy^iC) an d the ground state is estimated as 

\m'),ss^(0)\ < g^r-^mOWMCM- (87) 

Now we estimate ||Po^z(C)||- Since {^(Cj)}f=o * s a complete basis of the ground states, 
we can represent Pq^KC) by 

WHg^G), (88) 
where Cj is a complex coefficient. Thus, we have 

llwKC)ll 2 =Eii^ir(*'(0, *((,-)) = E C : 



3=0 



< ^ =i (*fD/ 3jr ^^| C7 .| j (89) 
where we have used and (|87jl. To evaluate ^ . |C 3 -|, we consider 

mCi),PoMQ) _£- r (JiQilM r i Vr (*(Ci),*(G)) rqm 

11^)11 ~ fe il*(0)llll*(C fc )ll J £j *||*(0)IHI*(Cfc)||" 1 ^ 

Then, we have 

1 jl \MCM f£ l kl \MQ\\\M( k )\\ 

< ^ l ( d+ m R -L^ + r -l*-i J2 \c k \, (91) 

where we have used (jHHj) and (|8*7jl. If we define |C m | = max{|Cfc|}^ , then we obtain 

G„ l(d+l)/2 R -L d - 1 

\C k \ < \C m \ < ^ i _ LdR _ Ld -, , 02) 
for any A; = 0, 1, • • • , L d . Thus, we obtain 

iip ^ mil 2 (^ + 1)^- 2 ^ 1 f93 , 

l|P0^(CJ|| <~Q- j _ L d R -L*-i • ^ 

from flU) and If we define L x by 1 - L x d Rr L ^ x = 1/2 and P 2 by P 2 = 2G 3 2 /G h 

then we obtain 

||P ^(C)I| 2 < F 2 (L d + l)l d+l R- 2Ld -\ (94) 
By definition of the normalized state (|66|). this inequality completes the proof of Theorem 
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4.4 Proof of Corollary 14.21 

Here, we prove Corollary 14.21 We define a normalized state ty± by 

*x := (1 - P °^ (1) . (95) 

||(i-p„)*l(i)|| 

An upper bound of (4f±,H4f±) gives an upper bound on the lowest excitation energy in 
a finite system, since ^j_(C) is orthogonal to all of the ground state. Since if^i(l) = 
H{1 — P )^i(l), the only remaining task is to estimate ||(1 — P )^l(1)||- From (JUJ, we 
obtain 

C 2 ( T d _|_ 1 "I T d + l G~2L d - 1 

||(1 - Po)*x(0)|| 2 =||*x(l)ir - II Wl(1)|| 2 < 1 - ^ [ iJ LdR L^ • ( 96 ) 
If we define L 2 by 

G, 2 (L 2 d + l)L 2 ^R-^ 1 

Gi 1 - Vit!-^" 1 2' 1 j 

and set P3 = 2i*i, then we obtain an upper bound 

(ty ± ,H^ ± ) < F 3 UL-\ (98) 

This gives an upper bound on the lowest excitation energy in a finite system, and so 
completes the proof of Corollary 14.21 I 



5 Existence of the Spin- Wave Gap 

In this section, we consider our model under the periodic boundary condition. The prop- 
erties of ground states and low energy excitations are very different from a system with 
the open boundary. We find only two ground states: the all-spin-up state and the all-spin- 
down state. We show that a one-magnon spin-wave excitation has an energy gap as in 
the XXZ model. The proof is based on Tasaki's argument for the SU(2) invariant model 
[TS] . He proved that the one-magnon spin- wave excitation in the Tasaki model has the 
same dispersion relation as that in the ferromagnetic Heisenberg model. These spin-wave 
excitations in both models have no energy gap, since they are the Goldstone mode above 
the ground states which spontaneously break the SU(2) spin rotation symmetry. On the 
contrary, an energy gap is generated by the anisotropy in our model as in the XXZ model. 
Here, we show only a brief sketch of the proof. 

5.1 Ground States and Spin- Wave Excitation 

First we obtain ground states. We have already found the representation of a ground 
state (|35p with the condition (|3bj) in section El The periodic boundary condition allows 
no configuration which satisfies the condition (|36|) except in the two cases: a x =| for all 
x G A Q or a x =[ for all x G A D . Thus, we conclude that all ground states in the periodic 
system are only two fully polarized states $| and 

Next we consider the one-magnon spin-wave excitation. For the spin-wave in our 
electron model, we consider properties of the spin-wave state in quantum spin models. 
The one-magnon spin- wave state with a wave- number k G /C: <3>sw(^) satisfies 

T x $sw(k) = e- ikx <S>sw(k) (99) 
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and 



'tot^SW — {Sxaax ~ J-J^SW, 

where x G A D . The translation operator T x is defined by 



T c T _1 

J -x^y,a ± x 



"x+y,o an d T x Cy a T x c x+y,a- 



K, is the space of wave-number vectors 

27m 



K :-- 



n e z a n 



L-l L-l 



Then, the one-magnon spin-wave state is in the following Hilbert space "H^ 



n k -.= \^en 



-ik-x 



* and S£}V = if I A 



1W 



(100) 
(101) 

(102) 
(103) 



We define the one-magnon spin-wave state with wave-number k by the lowest energy state 
in Tik. Let E$w(k) be the energy of one-magnon spin- wave state with wave- number k. 
We can prove the following theorem. 

Theorem 5.1 (Spin- Wave Gap) Suppose the d- dimensional Hubbard model defined by the 
Hamiltonian MoT) . There exist positive constants to, Uq, Xq, C < oo which are independent 
of system volume such that 



uanE sw (k)>^ 



d(\q\ + \q\- 1 -2) C 
2 A 



(104) 



for t > t ,U > Uq and X > X c 



This theorem shows that one-magnon spin-wave excitation has a finite gap for suffi- 
ciently large A in a periodic system. 



5.2 Sketch of Proof 



Theorem 15.11 is proved by the same approach given in ^H] • Here, we show only a sketch 
of the proof. 

First we introduce a localized electron operator a x a defined by 



V=£* 



V,a y,tr' 



(105) 



where ipy^l is defined by 



^ w := 



5. 



-p(tr)/4 



\ -x — e(j'),y "I - \ ".r+( 'J'./y 



+ 



eO'). 



A 



"x-eO\y + 



( 9 P(")/4)' 

A 



5 



x+e(i\y 



if x G A G 



if x G Aj 



(106) 



The set {a x ^vatJseA is a basis in the space of single electron states. We define the dual 
operator b x>a which satisfies 



{ b x,a, aJ >T } = 5 x ,ySa,r, and {6^^, 6 V)T } = = {a ] xa , aJ jT }. 



(107) 
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We represent 6 X (J in terms of the original electron operator by 

yeA 

Eqs. (Unnj), dTUTJ) and (fTUHl) mean 

E^S)*^ = ^ and E(®)*V'fc ) = (109) 

The original electron operators can be written in terms of a { xa and b X}(T , 

4,. = E(^)* a L and '•<•-• E'-'v- ( n °) 

The Hilbert space with |A D | electrons is also spanned by the basis 



A,B c A with \A\ + \B\ = |A | ^ , (111) 



because cj. CT can be written in terms of a\ a . 

We represent the interaction Hamiltonian in terms of this basis 



int 

xi , 



E f^E^W)^^) <, T <A,A 4 ,T- (H2) 

i,X2,a;3,a;4£A \ io£A / 



It is convenient to introduce a new hopping Hamiltonian H^ op defined by 

H hop :=t\ 2 E<- 6 ^ (ll 3 ) 

<j=t,J, xgA' 

for the estimation of a lower bound of spin-wave excitation, i^hop satisfies 

£y 1 1 (O ifxGA , 

Since i?hop<4 o-$vac = for x G A Q and t\ 2 is lowest energy eigenvalue of a single electron 
state which is orthogonal to the zero energy states, then we have i^hop < -f^hop- 

First, we define a basis of TCk- To define a convenient basis of Ttk, we define a state 
^^A{k) for /i = 0, 1, • • • , d and for a set i C A with |A| = |A„| — 1 by 

*m(*0 := E e ik ' WT A^,i ( II<T ) $ -c, (H5) 

weA \veA / 

where e M = o = (0, 0, • • • , 0) for fi = and e M = e® for fi = j (j = 1, 2, ■ ■ ■ , d). This 
state satisfies both properties and (jlUUJI . We define another state Q(k) by 

^ = ^HA E e *'™ T » a J At*t « *oa,\W, (US) 
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which is an approximation of the spin- wave state. We will choose a constant a(k) in the 
proof. We define the following basis of TCk by 

B k := {Q(k)} U | n = 0, 1, • • • , d, AC A 

with \A\ = \A a \ -1 and (fx, A) ^ (0,A o \{o})}. (117) 

We define H by H hop + H int and matrix elements h[$, between $, \& G Bk by the unique 
expansion 

M = h[W,$\V. (118) 

And we define D[$\ by 

D[$] := $]] - IM $ .*]I- ( 119 ) 

*eB fc \{*} 

Now, we can prove the following lemmas. 

Lemma 5.2 Lei E (k) be the lowest energy eigenvalue of H in the Hilbert space Tik- 
Then, we have 

E (k) > min DM. (120) 

<J>GB fc 



Lemma 5.3 There exist positive constants t , Uq, \q, C < oo independent of system vol- 
ume such that 



2U 



mm D[<^>] = D[Q(k)] > 



3=1 



c 



:i2r 



for t > to, X > \q and U >Uq. 



We find the proof of Lemma f5. 21 in subsection 6.1 of ref. |18j . Lemma f5. 31 is obtained 
by a direct evaluation. 

Now, we can prove Theorem 15.11 Since H < H, then we have E Q (k) < Egw(k). Thus 
we find 



2U 

AE > min En(k) > min -— 
fce/c fee/c A 4 



d(|g| + kl- 1 ) _ ^ eo S ( 2fc . e 0") + 6>) - g 

3=1 



:i22) 



from Lemma [5.21 and l5"31 This concludes Theorem 15.11 I 



6 Summary 

In this paper, we construct a set of exact ground state with a ferromagnetic domain wall 
structure and a spiral structure in a deformed flat-band Hubbard model under an open 
boundary condition. We have studied excited states above the domain wall ground state. 
There exists a gapless excitation above the domain wall ground state in dimensions higher 
than one. This excited state is constructed by acting a local operator near the domain wall 
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on the ground state. We study this model also under the periodic boundary condition. 
In this case a ground state becomes the all-spin-up or -down state. We have shown the 
energy gap of the spin-wave excitation above the all-spin-up ground state. These proper- 
ties of the excitations above the ground states are similar to those in the XXZ model. 
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